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Abstract 

Hopf algebraic structures will replace groups and group representations 
as the leading paradigm in forthcoming times. K -theory, co-homology, 
entanglement, statistics, representation categories, quantized or twisted 
structures as well as more geometric topics of invariant theory, e.g., the 
Grafimann-Cayley bracket algebra are all covered by the Hopf algebraic 
framework. The new branch of experimental mathematics allows one to 
easily enter these fields through direct calculations using symbolic ma- 
nipulation and computer algebra system (CAS). We discuss problems 
which were solved when building the BIGEBRA package for Maple and 
CLIFFORD ^ to handle tensor products, Grafimann and Clifford algebras, 
coalgebras and Hopf algebras. Recent results showing the usefulness of 
CAS for investigating new and involved mathematics provide us with ex- 
amples. An outlook on further developments is given. 



1. Aim of the paper 

The first aim of this paper is to present the features of BIGEBRA, a Maple V 
rel. 5 package. BIGEBRA was built to deal with tensored Clifford and Grafimann 
algebras, and it relies on CLIFFORD (Ablamowicz and Fauser, 1999-2002, 1996- 
2002). While the emphasis here is more on the package, we nevertheless provide 
novel results in the last section on quantum Yang-Baxter equations derived from 
a Clifford bi-convolution. CLIFFORD and BIGEBRA are meanwhile available for 
Maple 6, 7 and Maple 8. 

*Package homepage is located at http://math.tntech.edu/rafal/ while CLIFFORD 
and BIGEBRA are available for Maple V, Maple 6, 7, and 8. Maple ^ is available from 
http : //www .maplesof t . com. 

^The CLIFFORD package is described in a companion paper by Ablamowicz and Fauser 
(2003) 
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Secondly, we will not pass an opportunity to raise the flag for the subject of 
experimental mathematics. We strongly believe that experimental mathematics 
based on algorithmic approach has already been been changing research and 
teaching of mathematics and theoretical physics (Fauser, 2002). The main as- 
pects of experimental mathematics, in our view, are as follows: 

Computations 

Various results have been achieved by brute force calculations. CAS simpli- 
fies difficult computations with non-commutative algebras, and it allows one to 
tackle problems impossible for hand calculations. It makes fewer errors and its 
results can be validated. A CAS cannot and does not replace knowledge of the 
mathematics behind the problem and it requires that all concepts be well defined 
and sound. 

Checking Assertions 

Having a CAS at hand, one can easily check one's own assertions. A single 
counter-example uncovers one's misconceptions and leads to a more sound un- 
derstanding of the subject. Indeed, when checking a theorem on examples, one 
can occasionally find that it is not valid or that it has not been correctly for- 
mulated. This in turn helps with finding errors in the computer code (Never 
trust any result generated by a computer!) and/or with coming up with a proper 
formulation of such ill-stated theorem. 

Developing New Mathematics 

It has proved to be necessary to develop new mathematics for solving problems 
in physics. For example, we have introduced the Hopf Gebra as a weakened form 
of a Hopf Algebra. Using a CAS this was done by first exploring the set of prin- 
ciples which yielded correct physical results and then by fixing the mathematics 
of the new structure. 

Teaching 

Having a well developed CAS, it can be useful in teaching! Students can check 
their own prejudices by exemplifying them directly on a computer. This requires 
a stable code which does not allow for illegal input, etc. CLIFFORD' s code is stable 
while BIGEBRA still needs knowledge of the user because to gain speed, type 
and input checking is hardly done since it is time consuming. Hence, CLIFFORD 
(and with restrictions BIGEBRA) can be effectively used for a fast access to the 
mathematical field of GraBmann and Clifford algebras and cogebras in teaching. 
Students can make their own experiments and develop -given a CAS which is 
stable under silly input- a sound understanding of the mathematical structures 
in question. Teaching mathematics or physics using a CAS allows one to abstract 
from technical details on the first approach. However, if the goal is that students 
develop an understanding of the topic, they should be encouraged to recalculate 
by hand and to recode themselves. Otherwise no deep understanding will ensue. 
A good example how this can be achieved is Wright (2002). 
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Experimental Mathematics 

Having the opportunity to deal with a CAS opens the field of experimental 
mathematics. This includes partly the other topics given above but it should not 
be underestimated due to its own dynamics. Exploring mathematics by doing 
particular experiments and by strengthening or weakenning once's own assump- 
tions is of extreme value as it allows one to enter a new mathematical field 
quickly and in a secure way. 

Algorithmic Understanding 

We share a growing belief (Gruel and Pfister, 2002) that computer usage 
not only allows one to successfully complete non-trivial computations but also, 
through the developments of algorithms, contributes to a better understanding 
of the problem. The so called algorithmic approach leads to achieving a compu- 
tational proficiency that in turn leads to a deeper mastery of the subject on the 
theoretical level. 

2. The BIGEBRA package 

We concentrate in this article on BIGEBRA and its features which generally go 
beyond those of CLIFFORD and Maple. We assume that the reader is at somewhat 
familiar with the syntax and features of Maple, see, e.g., Wright (2002), while 
CLIFFORD is described in (Ablamowicz and Fauser, 2003). 

The package loads in a non verbose (silent) mode with the following commands 
while the linear algebra package linalg is loaded for convenience only. BIGEBRA 
automatically loads CLIFFORD which is needed for its functionality: 

i > restart : _CLIENV LSILENT] :=true : with (linalg) :with(Bigebra) ; 

[&cco, Szgco, &zgco-d, $zgpl-CO, &zmap, &u, EV, VERSION , bracket, cco_monom, contract, 
do_compose, dropj,, eps, gantipode, gco-d-monom, gcojmonom, gcojunit, gpl-cojmonom, 
gpLco-monom2 , gswitch, init, linop, Unop2, lists2mat, Iists2mat2, make.BIdd, mapop, 
mapop2, meet, op2mat, op2mat2 , pairing, peek, poke, remove_eq, switch, tcollect, tsolvel , 

type/ tens orbasmonom, type/ tens ormonom, type / tensorpolynom] 

The output is a list with available functions; some of them are for internal use 
only. CLIFFORD and BIGEBRA come with an extensive online help page system 
which is included in the Maple online help and can be searched. It contains 
not only the syntax of the procedures but even a good deal of unpublished 
mathematics. It is meant to present to the user the mathematical concepts at 
hand. Help topics may be reached in Maple fashion with the command ?Bigebra 
followed by the enter key. 

2.1. Tensor product 

Given the functionality of CLIFFORD to compute with Clifford algebras, we need 
two more key features to enter the realm of bi- and Hopf algebras. The first 
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one is the tensor product and some functions to manipulate it with while the 
second is the coproduct. The latter, however, follows naturally from the algebra 
structure on the linear dual space (Milnor and Moore, 1965). 

Maple comes with a define facility which allows to introduce ampersand 
operators &<name> that are associative (flat), commutative (orderless), linear or 
multilinear. This facility unfortunately has two drawbacks: 

• It cannot deal with user-defined scalars (ring elements). This is mathe- 
matically insensitive since any definition of a linear or multilinear function 
needs linearity with respect to "scalars" that one wants to compute with. 
We want to consider, e.g., tensor products (multilinear associative) over 
arbitrary rings, e.g., the integers, polynomial rings etc. 

• It produces code which gives wrong output. 

Remembering that GraBmann basis multivectors are denoted as Id, el, e2, ei, 
elwe2, elwe3, eiwej , . . ., let us explore properties of the associative (flat) mul- 
tilinear operator &r (comments in the code are separated by #): 

2 > restart : # unload BIGEBRA ! 

3 >def ineC&r' , flat, multilinear) : # produce the operator '&r' 

4 >out[l] :=&r(2*el,2*e2+4*a*e3) : # a is not known to be a ring element 

5 > out [2] : =el+e2 &r e3, &r(el+e2,e3) :# infix form is dangerous 

6 > out [1] ; out [2] ; # we delete output for typographical reasons 

4(el &r e2) + 8 (el kraeS) 
el + (e2 kr e3), [el kr e3) + (e2 kr e3) 

Notice that Maple does not know how to deal with a so, effectively, the mul- 
tilinearity is established over the integers only. If one wants to compute over a 
polynomial ring or, more generally, with K-modules, this function is useless. 

7 > constants : =constants, a: 

s >out[l] :=&r(2*e2,2*el+4*a*el) : # a is known to be a constant now 
9 >out [2] :=&r(true*el,false*e2) : # true and false are 'constants' 

10 > constants ; out [1] ; out [2] ; 

false, 7, oo, true, Catalan, FAIL, it, a 
A(e2 &r el) + 8 a (e2 kr el) 
true false (el kr e2) 

The last line does not make sense at all. Furthermore we find these peculiarities: 

11 >out [1] :=&r(-el,-e2) : # -> &r(el,e2) 

12 >out [2] :=&r(-el) : # -> -l*&r(el) 

13 >out[3] :=&r(2.5*el) : # -> 2.5*&r(el) 

14 >out[4] :=&r(&r(-el)) : # -> &r(el) 
is > out [1] , out [2] , out [3] , out [4] ; 
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(-el)kv(-e2), -{lkvel), 2.5(l&rei), -&r(l, 1, el) 

While the first expression is not even simplified, the last three results are just 
plain wrong. This error was detected in Maple 5 and is partially still present up 
to Maple 7. Hence, BIGEBRA includes a patched define facility. 

In CLIFFORD and BIGEBRA almost everything can be treated as scalars, hence 
one can compute over function spaces etc. The generic tensor product of BIGEBRA 
is given by &t. 

16 ># reload BIGEBRA, defined will be patched 

17 > restart : _CLIENV LSILENT] :=true:with(linalg) :with(Bigebra) : 

18 > out [1] : =el+e2 &t e3, &t (el+e2 , e3) : # infix form is still dangerous! 

19 > out [2] : =&t (&t (-el) ) : # fixed 

20 > out [3] : =&t (2 . 5*el) : 

21 >out [4] :=&t(a*el+2*e2,sin(phi)*e3) : 

22 > out [1] ; out [2] , out [3] ; out [4] ; 

el + (e2 &t e3), {el &t e3) + (e2 &t e3) 

-kt(el), 2.5&t(el) 
a sin(» (el &t e3) + 2 sin(^>) (e2 &t e3) 

The patched define which ships with BIGEBRA can handle tensor products over 
quite general rings. We give a final example where we define a tensor product 
&r over the polynomial ring Z[[x]] with integer coefficients. 

23 >def ineC&r' , flat , multilinear) : 

24 >unprotect ( 'type/cliscalar ' ) : # we want to change this type 

25 > 'type/cliscalar 1 :=proc(expr) type(expr,polynom(integer,x)) end: 

26 > out [1] : =type ( 2*x+3*x~5+5 , cliscalar) : # true, since in Z [ [x] ] 

27 > out [2] : =type( x/2 .cliscalar): # false, fraction 

28 > out [3] : =type ( 2*y~2+3*y-l , cliscalar) : # false, wrong indeterminant 

29 > out [1] , out [2] , out [3] ; 

true, false, false 

Now we are ready to use the new tensor product &r over Z[[a;]] : 

30 >out [1] :=&r(2.5*x,4*y) : 

31 >out[2] :=&r(3*x~2-x+5,x-x~4) : 

32 >out[3] :=&r(3*y*x~2-x-5,x*y-x~4) : 

33 > out [1] ; out [2] ; out [3] ; 

4x(2.5&ry) 

3x 3 (l&rl) - 3x 6 (l&rl) - x 2 (l&rl) + x 5 (1 &r 1) + 5 x (1 &r 1) - 5x 4 (l&rl) 
3 x 3 (y &r y) - 3 x e (y &r 1) - x 2 (1 &r y) + x 5 (1 &r 1) - 5 x (1 &r y) + 5 x 4 (1 &r 1) 

Observe, that real numbers like 2.5 or the variable y are still not treated as 
scalars but integers and x are. 
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2.2. Basic tools 

Having defined a tensor product, we need to have tools for operating on such 
structures. Let us fix notation as follows: A single term having no prefactors 
&t (bl , . . . ,bn) , where b\, . . . , b n are GraBmann basis multivectors, will be called 
a tensor basis monom or a word that is composed from letters of the GraBmann 
multivector alphabet. A tensor monom is a tensor basis monom with a "scalar" 
prefactor where "scalar" means a ring element, a function, a polynomial, or just 
a number. A tensor polynom, or just a tensor, is a sum of tensor monoms. It is 
the multilinearity that guarantees that every tensor can be written as a linear 
combination of some tensor basis monoms (Hilbert basis theorem). The i-th 
place in a list of arguments of a tensor will be called the i-th slot. 
Acting on a tensor means: 

• Removing or grafting of terms from or into a tensor 

• Rearranging the tensor 

• Acting with operators on n -slots producing m -slots 

The first set of operations is given by peek and poke. They work as expected: 
peek takes as an argument a number of the slot which it intends to remove 
and returns a sequence of lists of the removed terms and the remaining tensor. 
Procedure poke needs three arguments: a tensor, a GraBmann multivector which 
will be put into the i-th. place, and the slot number i. 

34 >f :=x-> 'peek' (x,2) :x:=&t(el,a*e2+b*e3,e3) :f (x)=eval(f (x)) ; 

peek(a&t(ei, e2, e3) + bkt{el 1 e3, e3), 2) = ([ae2, el &te5], [be3, el &t eS]) 



35 >g:=x-> 'poke' (x,e2,3) : x : =&t (el , a*e2+b*e3) :g(x)=eval(g(x) ) ; 

pokc(a(ei &t e2) +b(el kte3), e2, 3) =akt(el, e2, e2) + bkt(el, e3, e2) 

The second group of operations consists of switches -also called crossings or 
braids- which allow to reorder the tensors. From a CAS point, it would be 
easy, say, to multiply the i-th and the j-th slots of a tensor and place the 
output into the &-th slot. However, mathematical reasons disallow such a brute 
method. If one demands that any reordering process can be generated from the 
reordering rules of generators, one has to assume that a crossing is a natural 
transformation in a functorial sense and obeys coherence, see Kelly and Laplaza 
(1980); Lyubashenko (1995b, a). In fact these two conditions imply the quantum 
Yang Baxter equation which is discussed below. The GraBmann Hopf algebraic 
graded crossing and the ordinary switch (swap) of two elements fulfill these 
properties. 

The switch procedure swaps two adjacent tensor slots in a general tensor 
polynomial, while the graded switch gswitch respects the grading of the factors. 
Both functions take as a second argument i, the number of the tensor slot to 
act on, that is, to switch the i-th and the (i + l)-st tensor entries. 
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36 >fl:=x->' switch' (x,2) :x:=&t(el,e2,e3,e4) :f l(x)=eval(f l(x)) ; 

switch(&t(ei, e2, e3, e4), 2) = kt(el, e3, e2, e4) 



37 >f2:=x->'gswitch' (x,2) :x:=&t(el,e2,e3,e4) : f 2(x)=eval(f 2(x) ) ; 



gswitch(&t(ei, e2, e3, e4), 2) = -kt(el , e3, e2, e4) 



38 >x:=&t(el,e2,e3we4,e5) :f2(x)=eval(f2(x)) ; 



gswitch(&t(ei , e2, e3we4, e5), 2) ~ kt(el , e3we4, e2, e5) 

Notice the different signs for switch and gswitch depending in the latter mor- 
phism on the grading of the GraBmann multivector elements. 

Furthermore, BIGEBRA allows to act with any operator on tensors. We dis- 
tinguish such operations by the number of input and output slots. An ordinary 
endomorphism is a 1 — > 1 map while, e.g., a product is a 2 — > 1 map. Functions 
which are currently available are mapop and mapop2 for 1 — > 1 and 2^2 
operators, &map for 2 — > 1 operators and contract for 2 — > operators. Let 
us first show a linear operator projl acting on a certain tensor slot: 

39 >proj 1 : =proc (x) vectorpart (x, 1) end: 

40 >L:=x->' mapop' (x, 2, projl) :x:=&t(Id,el+elwe2,e3) :L(x)=eval(L(x)) ; 

mapop(&t(7d, el, e3) + kt(Id, elwe2, e3), 2, projl) = kt(Id, el, e3) 

where vectorpart (x, n) projects on the n-vector part of x in the GraBmann 
basis. Let us define a general element a; in a GraBmann or Clifford algebra over 
a two dimensional vector space and a general operator Rop given by its matrix 
R in a co-contravariant canonical basis: 



41 > restart : _CLIENV LSILENT] :=true:with(linalg) :with(Bigebra) : 

42 >dim_V:=2: # set dimension to 2 

43 >bas :=cbasis(dim_V) : # generate a list of basis monoms 

44 > X : =add(x [i] *bas [i] , i=l . . 2~dim_V) : # general element 

45 >R:=matrix(2~dim_V,2~dim_V, (i, j)->r[i,j]) : # matrix with entries r[i,j] 

46 >Rop:=proc(x) linop(x,R) end: # the operator Rop 

47 >X;evalm(R) ; 

X\ Id + x 2 el + x 3 e2 + x 4 elwe2 
n,i ri )2 t-1,3 r 1A 

r-2,1 r 2 ,2 7*2,3 7*2,4 
7*3,1 7*3,2 7*3,3 7-3,4 
_ 7-4,1 7*4,2 7*4,3 7-4,4 

The action of Rop on GraBmann multivectors and tensors is computed as 



48 >out[l] :=Rop(Id) :out[2] :=Rop(X) : out [3] :=mapop(&t(el,Id,e2) ,2, Rop) : 

49 > out [1] ; out [2] ; out [3] ; 
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n, i Id + ri, i el + r 3 , l e2 + r 4i i elwe2, 
(ri, i xi + n, 2^2+ ri, 313 + ri, 4 2:4) W + (r 3 , 414 + r 3 , 1 ari + r 3 , 2 2:2 + 3 x 3 ) e2 

+ (r2,2X 2 +r 2 ,3X3 + r 2 ,AX4 +r 2 ,iXi) el + (r 4; 1 x\ + r^ 2 x 2 +r 4;3 a;3 +r 4;4 x 4 ) elwe2, 

n, 1 &t(eJ , id, eg) + r 2 , 1 kt(el , eJ , eg) + r 3 , 1 &t(eJ , e2, e2) + r 4 , 1 kt(el , elwe2, e2) 

Operators may be defined as any Maple procedures also. 

The next class of operations is the application of product maps, i.e., 2 — > 1, 
which have two adjacent entry slots and one output slot. Examples of such 
operations which are predefined in BIGEBRA are the Grafimann and the Clifford 
products, and the contractions. User defined functions may be applied also. They 
are mapped onto tensors via the femap function which takes a tensor, the slot 
number and the (product) function as input parameters: 

bo >trm:=&t(Id,el,e2,Id) : ~ 

51 >out[l] := < &map(&t(Id,el,e2,Id) ,2, wedge) ' , ' — > 1 ,&map(trm,2,wedge) : 

52 >out[2] : = '&map(&t(Id,el,e2,Id) , 2 , cmul [K] ) ' , 1 — > 1 , &map (trm , 2 , cmul [K] ) : 

53 >out [3] : = '&map(&t(Id,el,e2,Id) ,2,LC,K) ' , ' —> ', &map(trm,2,LC,K) : 

54 > out [1] ; out [2] ; out [3] ; 

kmap(kt(Id, el, e2, Id), 2, wedge), — >, kt(Id, elwe2, Id) 
kmap{kt(Id, el , e2, Id), 2, cmul[K]), — >, kt{Id, elwe2, Id) + K h2 kt(Id, Id, Id) 
kmap{kt{Id, el,e2,Id),2,LC,K), >, K^ 2 kt{Id, Id, Id) 

Note that the Clifford product cmul [K] and the left contraction LC(...,K) de- 
pend on an arbitrary bilinear form K which is then passed on as an additional 
parameter to the procedure &map. This mechanism can then be used to deal 
with different Clifford algebras in different tensor slots. 

As a last example of actions on tensors we examine the evaluation and other 
2^0 mappings with values in the base ring. A very important case is given 
by the evaluation which employs the action of dual elements w.r.t. the canon- 
ical basis on multivectors. In fact one would need a new kind of basis vectors, 
however, for technical reasons of the current version of BIGEBRA, it is the user 
who has to take responsibility for keeping track of the tensor slots which are to 
contain co(multi)vectors. 

55 >out [1] :=[EV(Id,Id) ,EV(el,e2) ,EV(el,el)] : # eval = Kronecker delta 
se > out [2] : = ' contract (&t ( Id , el , e l+e2+e3 , Id) , 2 , EV) ' , 

57 > ' — > ' ,contract(&t(Id,el,el+e2+e3,Id) ,2,EV) : 

58 > out [1] ; out [2] ; 



[1, 0, 1] 

contract(kt(Id, el, el + e2 + e3, Id), 2, EV), — > , Idktld 
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2.3. Grafimann Hopf algebra 

A GraBmann Hopf algebra is a vector space that is a GraBmann algebra and a 
GraBmann coalgebra fulfilling certain compatibility laws (Sweedler, 1969). Since 
the algebra side is well known, we give some explanation on the coalgebra side. 
A coproduct on an algebra can be defined as the categorial dual of a product via 
a linear form defined on the vector space. This law is called product coproduct 
duality (Milnor and Moore, 1965). Now, it is easy to check that covectors may 
form a covector GraBmann algebra V v , where V is the GraBmann exterior 
product on co-multivectors. Using duality via the evaluation map (. | .) we find 

(w\xAy) = (w(i) | y)(w(2) I x), 

(W V w' | x) = {W | X(2))(w' I (1) 

where we have used Sweedler's notation for the coproduct A(x) = Yl( x ) x (i)® x (z) 
(the summation symbol is usually dropped). The evaluation map is given by 
the function EV from the package. The coproduct is seen to be a split of the 
homogeneous multivectors into pairs where the sign of permutation is taken into 
account. We give one example with symbolic indices a and b : 



>&gco(eaweb) ; 



(Id kt eaweb) + (ea kt eb) - (eb kt ea) + (eaweb kt Id) 

The first compatibility law valid in a Hopf algebra is that the product is an 
coalgebra homomorphism and the coproduct is an algebra homomorphism, see 
Fauser (2002). Furthermore, an antipode exists, which in the GraBmann case 
turns out to be isomorphic to the grade involution. We check this in dimension 2: 



>dim_V:=2: 

> 'S~ < =matrix(2~dim_V,2"dim_V, (i , j ) ->EV(bas [i] ,gantipode(bas [j] ))) ; 

> < Grade_Inv < =matrix(2"dim_V,2~dim_V, (i, j)->EV(bas [i] ,gradeinv(bas [j]))) ; 







-1 

1 



Grade_Inv = 





-1 
1 



The above described tools allow us to perform all algebraic manipulations in a 
GraBmann Hopf algebra, also using symbolic indices. 



2.4. Clifford bi-convolution 

A Clifford bi-convolution is a space endowed with a Clifford algebra structure 
and a Clifford coalgebra structure. Both products are unital and associative 
but not all such convolutions posses an antipode. However if an antipode exists 
then one can prove that one deals then with a Clifford Hopf algebra (Oziewicz, 
2001a; Fauser, 2002). We are interested in Clifford algebras over a general bilinear 
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form, hence we can introduce the Clifford product via the Chevalley deformation 
(Chevalley, 1997). Let x be an element in V, the space of generators, and let u 
be a general element in V A {= /\ V). Then, one defines the action of x on U as 

7 Z o u = i x (u) + x A u 

where i x is the inner product or contraction. The action is then extended to a 
map o : V <g> V i— > V by demanding that 

i x (y) = B(x, y), i u /w{w) = i u (iv(w)) 

i x (u A v) = i x (u) A v + u A i x (v) 

with u = (—l) 9u u being the grade involution, e.g., Fauser (2002). It is a re- 
markable fact that this product can be directly defined via a deformation called 
cliff ordization by Rota and Stein (Rota and Stein, 1994). It may also be called 
a Drinfeld twist, see below. Using the undeformed GraBmann product A and 
coproduct A one can write 

u o v = B A (u (2) ,v (1) ) u (1) A v {2 ) 

Here B A is the scalar valued bilinear form tied to the inner product action via 
the counit e as B A (a,b) = e(i a (b)). We exemplify this in the BIGEBRA package 
as 



63 > restart :_CLIENV[_SILENT] :=true: with (Bigebra) : 

64 >unprotect (gamma) : 

65 > Gamma :=proc(U) local x; x:=op(procname) ; LC(x,U)+wedge(x,U) end: 

66 >out [1] :=gamma[el] (Id)=Gamma[el] (Id) : 

67 >out [2] :=gamma[el] (e2)=Gamma[el] (e2) : 

6s >out[3] :=gamma[el] (e2we3)=Gamma[el] (e2we3) : 
69 > out [1] , out [2] , out [3] ; 

7ei(/rf) = el, 7 e i (e2) = B lj2 Id + elwe2, j el (e2we3) = B x . 2 e3 - Bi j3 e2 + elwe2we3 

If we do the same using the cliffordization process, which we will make explicit 
by defining a function cliff as to avoid using the internal function cmul, we get 

to > cliff :=proc(x,y) clicollect(simplify(drop_t(&map(contract( 

71 > &t (fegco (x) ,&gco (y) ) , 2 , scalarpart@LC) , 1 , wedge) ) ) ) end: 

72 >out[l] :=el &C e2we3 = cliff (el ,e2we3): 

73 >out[2] :=elwe2 &C e2we3 = cliff (elwe2 , e2we3) : 

74 >out[l] ;out[2] ; 

el kC e2we3 = B h2 e3 - B h3 e2 + elwe2we3 
elwe2 &C e2we3 = B 2 , 2 elweS - B x . 2 e2we3 - B 2 . 3 elwe2 - (-B 2 . 2 B 1>3 + B 2 . 3 B x . 2 ) Id 

The most interesting structure which we can now access computationally with 
the BIGEBRA package is that of a Clifford bi-convolution. In any pair ( U, V) of 
structures having a product on V <E> V and a coproduct on U, which need not 
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to be a bialgebra or Hopf algebra at all, one can define a convolution product 
between morphisms f,g,... : U — > V, see Fauser (2002). We consider here the 
GraBmann Hopf algebra and endomorphisms, and define the star or convolution 
product of endomorphisms according to Sweedler, (Sweedler, 1969) as 

(f*g)(x) = mo (f ®g) oA(i) 

where o is the composition of maps. The element x can be dropped safely. 
Re-examining the definition of the antipode shows that the map is the convolu- 
tive inverse of the identity morphism Idy on V. In the following steps we will 
compute the antipodes of the GraBmann Hopf algebra and of a twisted GraB- 
mann Hopf algebra where only the product is deformed, and of the Clifford 
bi-convolution, where both structure maps, product and coproduct have been 
deformed. The last possibility is beyond current deformation theory, where only 
one structure map is deformed. We choose dimension 2 for the generating space 
which gives a 4 dimensional GraBmann resp. Clifford algebra. The output is sup- 
pressed in most steps since its very tedious and long. The defining equations for 
the three antipodes are stored in eq_gr, eq_cl and eq_bc. 

75 > dim_V : =2 : bas : =cbasis (dim_V) : 

76 >B : =matrix(dim_V,dim_V, (i, j)->b[i, j] ) : # scalar product 

77 >BI:=matrix(dim_V,dim_V, (i, j)->p[i,j]) : # coscalar product 

78 >make_BI_Id() : ## <== initialize the Clifford coproduct 

79 > S := matrix(2~dim_V,2~dim_V, (i, j)->s [i, j] ) : # antipode template 
so >Sop := proc(x) linop(x,S) end: # operator using S 
si >X := add(x [i] *bas [i] , i=l . . 2~dim_V) : # general element 

82 > eq_gr : =drop_t (&map(mapop(&gco(X) , 1 , Sop) , 1 , wedge )-gco_unit (&t (X) , 1) ) : 

83 >eq_cl : =drop_t (&map(mapop(&gco(X) , 1 , Sop) , 1 , cmul ) -gco_unit (&t (X) , 1) ) : 

84 >eq_bc :=drop_t(&map(mapop(&cco(X) ,1, Sop) ,1, cmul ) -gco_unit (&t (X) , 1) ) : 

In a second step we solve these equations using the BIGEBRA tangle solver 
tsolvel: 



85>sol_gr : =tsolvel (clicollect (eq_gr) , [seq(seq(s [i , j] , 

86 > i=l. .2~dim_V) , j=l. .2~dim_V)] , [seq(x[i] ,i=l. .2~dim_V)]) : 

87>sol_cl : =tsolvel (clicollect (eq_cl) , [seq(seq(s [i , j] , 

88 > i=l. .2~dim_V) , j=l. .2~dim_V)] , [seq(x[i] ,i=l. .2~dim_V)]) : 

89>sol_bc : =tsolvel (clicollect(eq_bc) , [seq(seq(s [i , j] , 

90 > i=l. .2~dim_V) ,j=l. .2~dim_V)] , [seq(x[i] ,i=l. .2"dim_V)]): 



Before we display these morphisms in a matrix form, we compute a normalization 
factor for the antipode of the Clifford bi-convolution to simplify the output, 
which will be given as Sb c = N * S' bc . Here the primed antipode is the actual 
one and the unprimed is the normalized one. 



91 >N:=linalg[det] (linalg [diag] (l$dim_V)- B &* BI) : 

92 >S_GR := subs(sol_gr [1] ,evalm(S)) : 

93 >S_CL := subs(sol_cl[l] ,evalm(S)) : 
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94>S_BC := map(simplify@expand,subs(sol_bc[l] ,N*evalm(S) ) ) : 

95 >S_GR=evalm(S_GR) , S_CL=evalm(S_CL) ,S_bc=map(simplify,evalm(S_BC)) ; 



S.GR = 



1 











-1 

1 



S.CL 



S.bc = 



1 - Pi, 2 h, 1 + P2, 1 b 2 , 1+Pl,2 bi, 2 - Pi, 1 2 




Pi, 2 -P2, 1 






bi,2 


-62, 


1 
















1 









1 







1 












b\,2 


— 




1 













-1 















1 



Proposition 2.1: fdim_V=2 y ) T/ie antipode Sb c of a Clifford bi- convolution is 
up to a factor identical with the antipode S A of the Grafimann Hopf algebra if 
and only if the cliff ordization is performed with a symmetric scalar product and 
a symmetric coscalar product (both derivable from a quadratic and a coquadratic 
form via polarization in char 7^ 2) . 

COROLLARY 2.1: The recursive formula for computing the antipode given by 
Milnor and Moore (Milnor and Moore, 1965) and currently used in the Connes- 
Kreimer renormalization procedure cannot be applied in general for Clifford bi- 
convolutions with antisymmetric part in the scalar and coscalar product 

S(x) = e(x) — x — S(x'^)x'^, S(Id) = Id 

where the primed coproduct is over proper cuts only, i.e. x',^ 7^ Id. 

Note however that bilinear forms having antisymmetric parts are involved in the 
process of Wick normal ordering in QFT (Fauser, 2001). A CAS is a valuable 
help here to find explicit examples of such structures. 

3. Deformation and quantum Yang-Baxter equation 

In this section we want to use BIGEBRA to make the relation between cliffordized 
products, coproducts and standard deformation theory explicit. We keep the 
definitions from the previous section, i.e. dim_V=2, the settings for the scalar 
product B, the coscalar product BI, and the general element X. For further ref- 
erence we need an explicit form of the scalar product extended to V and recall 
the matrix form of the GraBmann antipode, this time from the BIGEBRA builtin 
function gantipode. We compute the convolutive inverse of the bilinear form 
needed below: 



96 >BW :=matrix(2~dim_V,2~dim_V, (i, j)->EV(cmul(bas [i] ,bas[j]) ,Id)) : 

97 >S_Gr :=matrix(2~dim_V,2~dim_V, (i, j)->EV(bas[i] , gantipode (bas [j] ))) : 

98 >BS :=evalm(BW &* S_Gr) : 

99 >BW=evalm(BW) , S_gr=evalm(S_Gr) ,BS=evalm(BS) ; 
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BW 



1 

&1,1 &1,2 

&2.1 &2.2 





BS 







&2,1 &1,2 — &2,2 &1.1 



1 













-&1,2 





-&2,1 


-&2,2 














-1 
0-10 

1 



&2, 1 &1,2 — ^2,2 &1, 1 

Note that the convolutive inverse is just the scalar product w.r.t the negative 
bilinear form. 

Proposition 3.1: Not every endomorphism has as the convolutive inverse B s = 
Bo S. However, exponentially generated morphisms and, more generally, graded 
morphisms do. 

A quasi triangular structure is an element R <E V ®V which satisfies, among 
others, the following condition: 



sw o A C i(x) = (%) <g> i? (2 )) o A Gr (x) 



(*) 



see, e.g., Majid (1995). Note that our definition is somewhat different due to our 
reversed ordering of tensor products of dual elements. The task is now to define 
a general element R and try to compute all possible solutions of equation (*). 
Therefore we start to define R as 



100 >Y:=add(y[i]*bas[i] ,i=l. .2~dim_V) : 

101 >R:=proc(x,y) local bas,tr_tbl; option remember; 

102 > bas : =cbasis (dim_V) : 

103 > tr_tbl:=table([seq(op(Cliff 5 [extract] (bas [i] ) )=i , i=l . .2~dim_V)]) ; 

104 > R[tr_tbl [op (Cliff 5 [extract] (x))] ,tr_tbl [op (Cliff 5 [extract] (y))]] 

105 > end : 

106 >RR:=add (add (contract (&t (bas [i] ,bas[j]) , 1 ,R) *&t (bas [i] ,bas[j]) , 

107 > i=l. .2~dim_V) , j=l. .2~dim_V) : 



Now the actual computation starts by evaluating the l.h.s. and r.h.s. of (*) as 



108 > Leq: =gswitch(&cco (X) , 1) : 

109 >Req:=&map(&map(switch(switch(&t(RR,&gco(X)) ,2) ,1) ,3, wedge) ,1, wedge) : 

Note that in the right hand side of the equation Req, two switches were used to 
put R in-between the xu\ of the coproduct and are not part of the formula (*) 
while the graded switch in the l.h.s. equation Leq is generic. This time we show 
directly how to solve tangle equations: 



no > eq: =tcollect (Leq-Req) : 

in > vars : =seq(seq(R[i , j] , i=l . . 2~dim_V) , j=l . . 2~dim_V) : 

112 >T:=seq(seq(&t(bas[i] ,bas[j]) ,i=l. .2~dim_V) ,j=l. .2~dim_V) : 

113 > CO : =seq(x [i] , i=l . . 2~dim_V) : 

114 > sys : = : 
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115 >for t in T do 
lie >for co in CO do 

ii7 > sys:=op(sys) , coef f (coef f (eq,t) ,co) ; 
us > od : od : 

119 >sol:=solve(sys,vars) : 

120 >matR:=matrix(2~dim_V,2~dim_V, (i, j)->R[i, j] ) : 

121 > 'R'=subs(sol,evalm(matR)) ; 



R = 





-Pl.2 







-P2.1 
-P2.2 









P2, lPl,2 - P2,2Pl, 1 



We have put the solution once more in matrix form. Inspection of this result 
shows, that our R is the convolutive inverse of the coscalar product, and hence 
it is exponentially generated. In fact it is an axiom of a quasi triangular structure 
that a convolutive inverse exists. We can check our assertion explicitly via 



122 > eql : =gco_unit (gco_unit (&t (X, Y) , 2) , 1) : 

123 > eq2 : =simplify (drop_t (contract (contract (&t (&gco (X) , fegco (Y) ) 

124 > , 2,scalarpart@cmul [-K] ) , 1 , scalarpartOcmul [K] )) )*Id: 

125 > is(eql=eq2) ; 



true 



Corollary 3.1: The (co-) cliff ordizations with respect to (co-) bilinear forms 
—K and K are convolutive inverse w.r.t. the Grafimann Hopf convolution. 

We are now ready to check that our quasi triangular structure fulfills the 
quantum Yang Baxter equation. We do this first for the graded switch of the 
GraBmann Hopf convolution, which is trivially a braid. We compute 

126 >Z:=add(z[i]*bas[i] ,i=l. .2~dim_V) : eqO:=&t(X,Y,Z) : 

127 > eql : =gswitch(gswitch(gswitch(eqO ,1) ,2) ,1) : 

128 >eq2:=gswitch(gswitch(gswitch(eq0,2) ,1) ,2) : 

129 > is(eql=eq2) ; 

true 

where we have not shown the cumbersome terms but simply checked that both 
sides evaluate to the same result, which establishes the assertion. 

Since a quasi triangular structure is an exponentially generated endomor- 
phism, we have to check that it fulfills 

R12R23R12 = -^23-^12-^23- 

To check this assertion we calculate with BIGEBRA : 



130 >bw:=proc(x,y) local i,j; 

131 > add(add(BW[i, j] *EV(bas [i] ,y) *EV(bas [j] ,x),i=l. .2-dim_V) , j=l. .2-dim_V) 

132 > end : 
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133 >Bsw:=proc(x,i) 

134 > tcollect(gswitch(contract(&gco(&gco(x,i+l) ,i) ,i+l,bw) ,i)) 

135 >end: 

136 >eq3:=Bsw(Bsw(Bsw(eqO,l) ,2) ,1) : 

137 >eq4:=Bsw(Bsw(Bsw(eqO,2) ,1) ,2) : 

138 > is (eq3=eq4) ; 

true 

which proves our claim in dim_V=2. As a last demonstration, we check the re- 
maining properties of a quasi triangular structure to be valid for an exponentially 
generated endomap. 

R(S(a),b) = R S (a,b), R s (a,S(b)) = R(a,b), R(S(a), S(b)) = R(a,b). 

Hence we generate the endomap 'R' and apply to its arguments the GraBmann 
antipode: 

139 >R: = 'R' : 

140 > out [1] : ='R'=matrix(2~dim_V,2~dim_V, (i , j) ->scalarpart (cmul [R] (bas [i] ,bas [j] ) ) ) : 

141 >out [2] :='RS c =matrix(2~dim_V,2~dim_V, (i, j ) ->scalarpart (cmul [R] 

142 > (gantipode(bas [i] ) ,bas [j] ) ) ) : 

143 > out [1] ; out [2] ; 



# 



RS 



"10 

Ri.i Ri, 2 

#2,1 #2,2 

#2, 1 Rl, 2 — R2, 2 Rl, 1 

10 

-#1,1 -#1,2 

— R2, 1 — #2,2 
#2, 1 Rl, 2 ~ #2, 2 Rl, 



144 >out [3] := < RSS < =matrix(2~dim_V,2~dim_V, (i, j ) ->scalarpart (cmul [R] 

145 > (gantipode(bas [i] ) ,gantipode(bas [j] ) ) ) ) : 

146 >out [4] := < SR'=matrix(2"dim_V,2~dim_V, (i, j ) ->scalarpart (cmul [R] 

147 > (bas[i] ,gantipode(bas[j])))) : 

148 > out [3] ; out [4] ; 



RSS 



SR 





1 
















#1, 1 


#1,2 










#2,1 


#2,2 
















#2, 1 #1, 2 — #2, 2 #1, 1 


1 


















#1,1 


— #1,2 










#2, 1 


-#2.2 
















#2, 1 #1,2 - #2,2 #1, 1 



Comparing the matrix representations gives us the desired results when we recall 
that the convolutive inverse is given by the 'scalar product' -R on the generating 
space. 
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We can compute along the same lines as exemplified above the Yang Baxter 
matrix, a 16 by 16 matrix, which represents the action of R on V <8> V. Since 
the output of these computations is quite lengthy we will give only some further 
results and not the explicit matrix. The actual worksheet is available from the 
home page of the second author. 

Remark: The Yang Baxter matrix of our 2 dimensional example has determinant 
1 and is a function of all four parameters of the quasi triangular structure R. 

A further question is, if this structure is really quasi triangular or only trian- 
gular. Due to our reversed indexing in duals, we have to check for triangularity 
if the Yang Baxter matrix squares to unity. 

Remark: The Yang Baxter matrix of our 2 dimensional example is triangular if 
and only if the exponentially generated endomap is derived from a symmetric 
bilinear form on the generating space. This has deep implications for ordering 
process in quantum field theory, see Fauser (2001). 

We expect such assertions to be true in higher dimensions, though that needs 
an algebraic proof. However, having solutions in low dimensions allows in many 
cases to prove by induction a general hypothesis. This is the step from experi- 
mental mathematics to 'pure' mathematics. 

4. Conclusions 

The present paper was intended to show how a CAS allows to ask questions of 
research interest and to come up with low dimensional solutions. During this 
process we have worked out examples showing how to use the BIGEBRA package 
and how to attack more advanced problems. 

We hope that it became obvious from our presentation of the CLIFFORD 
(Ablamowicz and Fauser, 2003) and BIGEBRA packages that the ability to com- 
pute, via a CAS, allows to build a sound knowledge about mathematical prob- 
lems and that there is a need for experimental mathematics in education and 
research. 

This article is, due to restrictions in space, obviously not able to give a com- 
plete review of all abilities of the CLIFFORD and BIGEBRA packages. Therefore 
the interested reader is invited to check the package home page* for the doc- 
umentation. Both packages come with a built-in online help where, for every 
function, syntax, synopsis and examples, and sometimes a more advanced math- 
ematical background is provided. A printed ps or pdf version of approx. 500 
pages is available there also. For those who wish to have a look at and a feel for 
the packages, there is a possibility to access them online over the web§. 

turl: http : //math . tntech . edu/raf al/ 

Wl: http://kaluza.physik.uni-konstanz.de/~ fauser 
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There are tremendously many open problems coming with a bi-convolution 
of mathematical and physical type. It is not yet clear, under which conditions 
antipodes exist in higher dimensions. In dimension 3 a Clifford bi-convolution 
has in general no antipode, hence further relations between scalar and coscalar 
product have to be found. Only very little is known about the nature of the 
crossings derived from antipodal bi-convolution, but see Fauser and Oziewicz 
(2001). The axiomatics, in terms of morphims and category theory of Clifford bi- 
convolutions, is under consideration but not yet finished (Oziewicz, 1998, 2001b). 
Clifford bi-convolutions are deeply connected to the renormalization process of 
the quantum filed theory. However, there one has to go beyond the scheme 
presented in this paper, see Brouder (2002); Brouder et al. (2002). A CAS and 
especially CLIFFORD and BIGEBRA are ideally suited to explore this exciting 
field. 

Finally we want to emphasize that the present article contains Maple output 
(with only a minor TgX cosmetics) which was generated directly by processing 
a Maple worksheet. 

Acknowledgment: The second author, BF, acknowledges gratefully financial 
support from University of Konstanz, LS Prof. Heinz Dehnen. 
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